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We study here the evolution of a massless scalar field in a spacetime, developing from 
a regular initial spacelike surface. The Einstein equations and regularity and boundary 
, conditions governing the same are specified. Both homogeneous and inhomogeneous collapse 

^ ■ models are considered and we analyze when the occurrence of singularity will be simultaneous 

^ ■ or otherwise. In the inhomogeneous collapse case, we characterize a wide family of black 
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hole solutions arising in scalar field collapse. We also discuss the possibility of existence of 
classes of non-singular models where collapse could almost freeze if suitable conditions are 
satisfied. 
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CN ■ I. INTRODUCTION 

> 

in 

I The formation of spacetime singularities in gravitational collapse and formation of black holes 

is an issue of great importance in gravitation physics, which has been investigated in much detail 

!>■ 

' in Einstein's theory. The occurrence of singularities offers the regime where gravity is extreme, and 

oo ■ 

' where the quantum gravity effects would be important. As is known, dynamical evolution of matter 

- ^ , fields in a spacetime generically yields a singularity, provided reasonable physical conditions are 

^ • satisfied such as the causality, a suitable energy condition ensuring the positivity of energy density, 

and formation of trapped surfaces. 

The case of gravitational collapse of a massless scalar field is of particular interest in both 
collapse situations as well as cosmological scenarios. In cosmology, special importance is attached to 
the evolution of a scalar field, which has attracted a great deal of attention in past decades. This is 
because one would like to know the behaviour for fundamental matter fields towards understanding 
the transition from matter dominated regime to dark energy domination (see e.g. [l| and references 
therein). Scalar fields are of much interest in view of the inflationary scenarios that govern the early 



universe dynamics because such a field can act as an 'effective' cosmological constant in driving the 
inflation In gravitational collapse studies, the nature of singularity for massless scalar fields 
has been examined and a number of numerical and analytical works have been done in recent years 



on spherical collapse models [3[-[10[], from the perspective of the cosmic censorship hypothesis. 



Hi, 



The massless free scalar field has been studied for the static case also in some detail ( see e.g. 
1^ . [l^ . jl^ - [isl . However, the dynamical case is more important to understand because 

that may give insights into phenomena such as gravitational collapse and the cosmic censorship 
hypothesis, other than the early universe and cosmological considerations. 

In the present study, we first develop here a mathematical structure in a general manner to 
deal with the evolution of massless scalar fields in a spacetime. This would be applicable to 
problems either in cosmology or for gravitational collapse final states in a spherically symmetric 
spacetime. We then examine the classes of collapsing models where the singularity in future forms 
simultaneously as the collapse develops, even when the density of the field could be inhomogeneous. 
We study here both homogeneous as well as inhomogeneous collapse models. We then characterize 
here a wide class of black hole models forming as collapse final state. Apart from the spherical 
symmetry, we do not assume here any further constraints on the spacetime, such as the presence of 
a homothetic Killing vector, homogeneous or shearfree nature for the fluid, or such other conditions. 

We note here the classes of models where our considerations would apply. The analysis in 
this paper has been done using a comoving coordinate system. For massless scalar fields, such a 
coordinate system would break down when the gradient of the scalar field becomes null. So the 
conclusions based on the present analysis would be valid for those classes of models and solutions 
for which the gradient of the scalar field is timelike throughout the dynamical evolution of the field. 
We note that the homogeneous and isotropic Friedman-Roberson- Walker solutions with a massless 
scalar field as the matter content is such an example, where this condition is satisfied. Massless 
scalar field solutions with inhomogeneous perturbations around a homogeneous background would 
also be in this class. In fact, as we shall argue later, this class where the gradient of the scalar field 
always remains timelike, would include a large number of physically relevant situations. On the 
other hand, our analysis does not apply to the case where the gradient of the scalar field changes 
its sign, becoming null and then spacelike, from its original nature of being timelike. However, as 
we have argued here, the present class is of sufficient physical interest in its own right to carry out 
an analysis of the same as far as the dynamical evolutions of the scalar fields are concerned. A 
related important point that we note is, the analysis given here would also hold and describe the 
evolution of stiff fluids in a spacetime in a general manner. This is because, a massless scalar field 
with a timelike gradient, which is minimally coupled to gravity, has an exact correspondence with 
a stiff fluid minimally coupled to gravity. 
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The massless scalar field (j){x'^) on a spacetime {M,gab) is described by the Lagrangian, 

= -^(p-a4>;bg'''' ■ (1) 

The corresponding Euler-Lagrange equation is (p-abd"'^ = 0, and the energy-momentum tensor for 
the scalar field, as calculated from the above Lagrangian is, 

Tab = (t^;aCt>,b " \gab (j^ A^d'^) ■ (2) 



A massless scalar field is a Type / matter field [17|], i.e., it admits one timelike and three spacelike 
eigen vectors. At each point q G M, we can then express the tensor T""^ in terms of an orthonormal 
basis (Eq, El, E2, E3), where Eq is a timelike eigenvector with an eigenvalue p and E^ (a = 1,2,3) 
are three spacelike eigenvectors with eigenvalues pa • The eigenvalue p represents the energy density 
of the scalar field as measured by an observer whose world line at q has an unit tangent vector Eg 
and the eigenvalues pa represent the principal pressure in three spacelike directions E^. We now 
choose the spherically symmetric coordinates {t,r,9,(j)) along the eigenvectors (Eo,Eq), such that 
the reference frame is comoving. As discussed in [isj, the general spherically symmetric metric in 
comoving coordinates can be written as, 

where dil^ is the metric on a unit 2-sphere and we have used the two gauge freedoms of two 
variables, namely, t' = f{t,r) and r' = g{t,r), to make the gtr term in metric and the radial 
velocity of the matter field to vanish. We note that we still have two scaling freedoms of one 
variable in t and r. 

In general, we have <j) = '/'(^j'^)) but from equation ^ it is easily seen that in the comoving 
reference frame ([3]), we must have (p{t,r) = (p{t) or 4'{t,r) = (j){r), because the energy-momentum 
tensor is diagonal. As we would like to investigate here the dynamic behaviour of the scalar field, 
we consider the former option. In this comoving frame the components of the energy-momentum 
tensor are. 

Thus, we see that in the comoving frame the massless scalar field behaves like a stijf isentropic 
perfect fiuid with the equation of state 

p(t,r)=p(t,r) = e-2-(*.'-)0V2. (5) 
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We can easily see that for any real valued function (j){t), all energy conditions are satisfied by 
the matter field. We note that if we consider the field (p = (j){t), then the weak energy condition 
guarantees that (p,^ is either timelike or null always in general. However, in a comoving frame, the 
gradient of the scalar field remains timelike throughout the collapse, which is a property of the 
comoving frame by definition. 

In a physically reasonable collapse situation, the energy density of the matter field and the 
magnitude of the Ricci scalar are expected to increase with time. For a massless scalar field we 
have I (p,^(j)'^ 1= 2/3. In such cases, if collapse progresses from a regular spacelike hypersurface 
where the gradient of the scalar field is timelike, then the density is non-zero initially and it can 
only increase afterwards. This would mean that throughout the collapse evolution of the scalar 
field from the given initial regular surface, the gradient of the scalar field would always remain 
timelike. So in this case of a collapsing scalar field, we can use the comoving coordinate system 
without any loss of generality and without a concern on a possible breakdown of the coordinate 
system. 

It is important to note that we are not specifying here the class of initial data from which one 
can have collapse situations where the density would always increase. But given such a solution, 
our analysis would hold. Therefore, the comoving system can be used to study the collapse of a 
massless scalar field, where the density either increases or does not decrease with time. In this 
paper, by a collapsing scalar field, we would mean only such cases. 

The correspondence here with a stiff fluid can be seen from ([2]). The energy momentum tensor 
for the perfect fluid is 

Tab = {P + p)UaUb + P9ab, (6) 

where is the velocity vector. For stiff fluid, p = p. 

Since (f),n is timelike, (h,a4>'^ = — I 4>,u4''^ I- Defining n„ = — "^"^ ^ , the energy momentum tensor 
for the massless scalar field can be expressed as. Tab = (I 0'^ D^a^b + ^dabil 4>jfi 4'''^ I)- 

Denoting | \= p = p, this expression for the energy momentum tensor is the same as 

that for the stiff fluid. The unit velocity vector, = {(l){t))~^ In a comoving coordinate system, 
we can choose li'^ = (1,0,0,0), which for the massless scalar field would be equivalent to choosing 
(f){t) = 1. If does not diverge, such a choice is always possible. 
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II. EINSTEIN EQUATIONS, REGULARITY AND BOUNDARY CONDITIONS 

The dynamic evolution of the initial data, as specified on a spacelike surface of constant time 
is determined by the Einstein equations. For the metric ([3]), using the definitions 

G{t,r) = e-^^{R'f,H{t,r) = e~^\R)\ (7) 

and 

F = R{l-G + H), (8) 

the independent Einstein equations for the massless scalar field (in the units SirG = c = 1) are 
then given by, 

F' = h-'^'^'^R^R! , (9) 
F = -^e-2>2i^2^ , (10) 
dt (fl2e^->) = , (11) 

R'G - 2Ru'G = . (12) 

Here (') denotes the partial derivative with respect to the coordinate r and (') with respect to t. 
The function F = F{t, r) has an interpretation of the mass function for the collapsing cloud, and 
it gives the total mass in a shell of comoving radius r on any spacelike slice t = const. The energy 
conditions imply F >0. 

The function R(t, r) is the area radius of a shell labeled r at an epoch t. For the sake of 
definiteness let us consider the situation of a collapsing cloud, and we have i? < as we are 
considering the collapsing branch of the solutions. If R changes sign then that corresponds to 
a bounce or dispersal of the field during evolution. We use the scaling freedom for the radial 
coordinate r to write i? = r at the initial epoch t = ti, and in order to distinguish the regular 
center of the cloud at r = from the genuine spacetime singularity at the termination of collapse, 
where the area radius i? = in both cases, we introduce a function v{t, r) as defined by 

v{t,r)=R/r. (13) 

We then have, 

R{t, r) = rvit, r), viti,r) = 1, u(t,(r), r) = (14) 
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with -u < 0. The time t = ts{r) here corresponds to the shell- focusing singularity at R = 0, where 
the matter shell labeled a comoving radius of constant r collapses to a vanishing physical radius R 
on reaching the genuine spacetime singularity. 

We note that equation (jlip is the Klein-Gordon equation for the scalar field, which is a part 
here of the Einstein equations via the Bianchi identities. We can integrate this equation to get 

R^e^-^cj) = r^f{r), (15) 

where /(r) is an arbitrary function of integration. We can now eliminate the function (j){t) from 
equations Q and ([TO]) to get 

F' ^ F ^l r^f{r)G 
R' R 2 R'^R'2 ■ 

We now have four Einstein equations, namely ([8|), (jl2p and (jl6p . and four unknown functions of 

two variables, ^, R and F. Solution of these equations, subject to the initial data and energy 

conditions, would determine the time evolution of the system. 

Our purpose now is to construct the classes of solutions to the Einstein field equations, which 

give the dynamical scalar field collapse evolutions, given the initial data at an initial time t = ti. 

We define the suitably differentiable functions M{r,v) and A{r,v) as below 

M{r,v)^^^, (17) 

Mr,v),,^^. (18) 

We note that Einstein equations ([9|) and (jlOp . imply that F must behave necessarily as closer 
to the regular center r = of the cloud, in order to preserve the regularity of initial data and to 
preserve the finiteness of matter density at all regular epochs of evolution. Hence, as M{r,v) is a 
general, at least function, the equation (jl7p is not really any ansatz or a special choice, but a 
fully generic class of mass profiles for the collapsing cloud, consistent with and as allowed by the 
regularity conditions. 

Also, in order to be specific, we construct only classes of collapse evolutions which admit no 
shell-crossing singularities in the spacetime where R' = 0. We therefore consider only the genuine 
singularity at = where the physical radii of the collapsing shells vanish, and not the cases 
when nearby shells of matter may cross, giving rise to a density singularity which need not be 
gravitationally strong. Therefore the function A{r, v) is well-defined for all non-singular epochs. 
Now using equation (llSp in equation (I12|) we get, as a class of solutions of Einstein's equations 

^(r,?;) =6(r)e2'^^('^''') . (19) 
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Here b{r) is another arbitrary function of the shell radius r. The regularity condition on the velocity 
function v at the center of the cloud implies that the form of b{r) has to be, 

b{r) = l + r^bo{r). (20) 

A comparison with the Lemaitre-Tolman-Bondi dust collapse models j3] implies that we can 
interpret 60 (^) ^-s the energy distribution function for the collapsing shells. 

We emphasize that the functions M{r,v) and A{r,v) are not independent here in the case of 
scalar field collapse. Because, using (I19p in the second part of ()16p gives the required relation 
between them as 



2rA{r, v) =ln 



-2M{r, v).uV^{v + rv')^ 



(21) 



where q{r) = /^(r)[l + r-^6o(r)]. Now to determine the function 7V4(r, u), we use ()17p in the first 
part of (fT6]) to get the required first order equation 



3A^(r, v) + rM{r, v)^r + Q{r, v)M{r, v)^^ = . (22) 
where Q{r,v) = {2rv' + v). The above equation has a general solution of the form 

HX,Y) = 0, (23) 



where X{r,v, Ai{r,v)) and Y{r,v,Ai{r,v)) are the solutions of the system of equations, 

dM.{r,v) dr dv 



(24) 



'iM{r,v) r Q 

Amongst all the classes of solutions M{r,v) of above, those are to be considered which obey the 
energy conditions and required regularity conditions for the collapse. 

To specify the boundary conditions for the cloud, solving (|22p at r = 0, we get 

lim7W(r,?;) = ^ , (25) 

where mo is a constant and from ()2ip we see that the regularity at the center of the cloud requires 
6mo = q(O). Along the singularity curve u = the density diverges. To solve for the function 
v{t,r), we use the equation of motion ([8]), and defining a function h{r,v) as 

h{r,v) = [e^''^^'''''> -l]/r^, (26) 

we get, 

^v = -g(r,v) (27) 
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where, 

g{r,v) = e^(^'^) ^J vboe^"^^ + vh{r, v) + M (r, v) . (28) 



The negative sign in the right hand side of equation (j27[) corresponds to a collapse scenario where 
we have R < 0. Also in the (r, v) plane, the function u(r, v) is related to the function A(r, v) by 
the relation 

i/(r, v)^r + '^{f, v)^vv' = A{r, f ),t>(f + rv') (29) 

Given the functions A{r,v)^v and v' in terms of r and v, this is again a quasi-linear first order 
partial differential equation like (|22p . with a similar general solution. 

Now we have derived the relations between all the functions in the (r, v) plane. The solutions to 
this system describes the evolving scalar fields which we discuss, and then we consider the nature 
of initial data required to completely specify the collapsing model. 

To see the existence of solutions to the complete system evolving from an initial spacelike slice 
t = ti, consider the partial derivatives of the function v{t, r). We get v' from ([22]) 



and from ([27]) we get v. 

. _ G{r,v) 



(31) 

\/v 

These give expressions of v and v' in terms of v,r, M(r,v) and its derivatives. As we have 
already seen from (|21|) . we can write the function A{r,v) and hence u{r,v) in terms of A4{r,v) 
and it's derivatives, r,v and v' . Thus, to get a solution of v{t,r), which would in turn provide the 
complete dynamical solution to the collapsing system, we need the Pfaffian differential equation, 
v'dr + vdt — dv = 0, to be integrable. The integrability condition is given by. 



. / 

'if 



(32) 



This integrability condition gives the required second order equation for the mass function M{r, v). 
Any solution of the above equation would by default solve the quasilinear equation (|22|) . via (|3Up . 
Also, this would then uniquely solve for functions A{r,v) and hence z/(r, ?;) via ()2ip . The solution 
set of this integrability condition is non-empty as Ai{r,v) = uiq/v'^ solves the whole system with 
V = v{t), V = v{t) and A = A{r) to give a Fi? 14^ interior. 

For each solution M of (1321) . we would have a particular Pfaffian equation, which would be 
integrable. We know that given one integrating factor for a Pfaffian differential equation, we can 
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find infinity of them. So for any sucli M, if we have one such integrating factor of the corresponding 
Pfaffian equation, then there would be an infinite number of such integrating factors. In general, 
they would correspond to an infinite number of solutions of that Pfaffian equation. Among the 
solutions that exist, we choose only those which obey the required boundary conditions. 

We now consider the independent initial profiles required at the epoch t = ti or v = 1 to evolve 
the system via Einstein equations. From ([9]) we see that providing the function 7W(r, 1), which 
is the same as F{ti,r)/r^, would determine the function v{ti,r) (upto a multiplicative constant 
(j){ti)). This would then determine the function A{r,v) ^y\y=i using (jlSp or (j29p . Also from (j22p . we 
see that 7W(r, 1) specifies A4(r, f )^„|„=i, which in turn determines A[r,l) using (12ip and h{r,l). 
Hence we see that specifying the functions 7W(r, 1) (interpreted as the initial density profile) and 
— + f'^^oi''')] (which specifies the initial velocity profile) completely determine the 

evolution of the system. 

Assured of the existence of solutions, integrating equation (j27p with respect to f , we can write 
the solution of v{t,r) in an integral form, 

Note that the variable r is treated as a constant in the above equation. The above gives the time 
taken for a shell labeled r to reach a particular epoch v from the initial epoch v = 1. 

One could sketch the iterative process by which the initial Cauchy data evolves. As we noted, 
given the density and energy profile we can find all the other functions at the initial epoch v = 1. 
Then we use (j33p to find the functional form of v{t,r) for v = 1 — e where e is an infinitesimally 
small positive number. As we easily see, this functional form would be 

v{At,r) = 1- j{r)At, (34) 

where 

j(r) = e"''^'^^ ^ boe-^''-Mr) + ho{r)+Moir). (35) 

Here the subscript denotes the initial data and At denotes the infinitesimal evolution of the 
spacelike initial slice in (t, r) plane. We now calculate the function v' as 

v' =j'{r)At = j'{r){l-v)/j{r) (36) 

Using this in (j22p gives M{r,v) at v = 1 — e with respect to the given initial data Aioir) and 
boundary conditions. Using the form of Ai{r,v) and v' we then calculate the function A(r, 1 — e) 
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using ([2T]) . i/(r, 1 — e) using ([29|) . and finally /i(r, 1 — e). We then plug all these functions in ([33|) 
again for the next iteration, until we reach v = 0. 

Now we see that the time taken for a shell labeled r to reach the spacetime singularity at -R = 
(which is the singularity curve), is given as 



In a physically realistic gravitational collapse situation such as collapse of a massive matter cloud 
which continually collapses, one would focus only on those classes of solutions where ts(r) is finite 
and sufficiently regular. This means that the cloud collapses in a finite amount of time. 

III. COLLAPSE OF THE SCALAR FIELD TO A SIMULTANEOUS SINGULARITY 

We now consider the case when the singularity occurring in the spacetime as a result of collapse 
is simultaneous. This would be the case, for example, when the density is same at every point 
in space at any given time, i.e. the collapse is homogeneous. In this case, p = p{t). Since 
p = e~^'^(/>(t)^, therefore we have u = v{t) only. We rescale t to make e'^'^ = 1. The spacetime 
singularity occurs when the physical radius goes to zero, i.e. v = 0. In the case of the density 
being homogeneous, this implies that the singularity curve ts{r) is independent of r. 

The expression for ts{r) in general is, 



Here the r-independence of ts{r) then implies that the integrand on the right hand side is a function 
of V only. The time coordinate can be expressed in general as 



where /ii(r) is an arbitrary function of r. As the integral is a function of v only, so the initial 
condition which is v = 1 at the time t = ti, implies that hi is a constant. Therefore from (j39p . it 
is seen that 




(37) 




(38) 




(39) 



v' = 0. 



(40) 



^From the Einstein equations, we get the relation 



e-2>(t) 



2 



2M. 



(41) 



V 



,2 
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Since the left hand side of the equation is a function of t only, the right side must be a function of 
v{t) only. This implies M = M{v) only. So ()22|) gives 

M = ^ (42) 

^Prom (fTSl) we get, 

A,^ = (43) 

This implies A = A{r). This is consistent with the form of G we get from the other equation, 

2v^M,, {v + rv'f 
G - (44) 

Putting M in the last equation, we get G = j?^- Since ts{r) ^ 0, the integrand in its expression 
must be finite at r = 0. This implies, 

where /i(0) is finite. In this case, since ts(r) is a constant, /i(r) is a constant also. So we can write 

f{r) = (46) 
where c is a constant. So in this case we have e'^^ = j^^- ^From (j39p we get, 



-dv (47) 



ds'^ = dt^ - v{tf{^^—^ + r^dO^) (48) 



The metric then is given as, 

dr 
1 + cr 

In this case, ts{r) is finite and (j){t) blows up at the singularity. This can be seen from (j4ip 
as it reduces to (j){t) = We also note that c = 0,±1, gives all the possible solutions in 

this case, which correspond to the Priedmann-Robertson- Walker models. When c = 0, we have 
v{t) = [Ki — Sy^mot)^, where Ki is some positive constant. 

We note here that we considered above the case when the density p = p{t), which leads to a 
simultaneous singularity, and the PRW class of models. However, even when the density is not 
homogeneous, a simultaneous singularity can result as we shall discuss below. 



IV. COLLAPSE OF A SCALAR FIELD OF INHOMOGENEOUS DENSITY 

In this case, the singularity that results as the collapse endstate need not be simultaneous in 
general, and in general we have p = p{r,t). Before proceeding further, we discuss the regularity 
conditions that are required for a physically reasonable collapse of the scalar field. 
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A. Regularity conditions 



First, we note that at the center when r goes to zero, we must have 



lim(rt>') = 



(49) 



Because, when this condition is violated, v becomes divergent as r goes to zero at the center. 
For example, let linir^Q^rv') = C2, where C2 has some non-zero value. Then liuir-^QV goes as 
mod Inr. If \iinr—,Q(rv') is divergent in r, then the divergence of v will only become more severe. 
But divergence of v means that R' becomes arbitrarily large for the comoving shells near the center. 
This is clearly unphysical and is to be ruled out. 

The second condition comes from one of the Einstein equations. In the limit of going to the 
center limr — > 0, using (j49]) in ([22]) . we get 



Here limr^o{rM,r) can behave in three possible ways, which are, i) limr^Q{rM,r) = 0, (ii) 
lim,r^o('"-^)r ) has some finite non-zero value, and (iii) \imr^o{rM,r ) can be divergent. 

Now we point out that only the first option is possible for a physically reasonable collapse 
model. In the case of (ii) occurring, we have limr^o M = co(/nr) + h2{v). Now we can argue that 
h2{v) is not divergent. Otherwise, M,^ will also contain a logarithmic divergence. But for v ^ 0, 
this implies a divergence of density, which is not physical. So we cannot consider that case. For 
the other case, one can see from (I38p . that the integrand becomes divergent when r goes to zero. 
Therefore ts(0) = 0, which means that the singularity is present even at the initial epoch. So (ii) 
is ruled out. Similarly, one can also rule out the possibility (iii), proceeding in the same way as for 
the case (ii). In this case also M has r-dependent divergence, which makes ts{0) = 0. Therefore, 
(i) is the only possibility that is allowed. 

In such a case then, (]50p reduces to 3M + vM,^= in the limr ^ 0. This implies 
lim^^o ^(^'j ^^) = rrio/v^, for -u / 0. 



3M + rM,r +vM,y = 



(50) 



B. Black hole formation in scalar field collapse 



In the following, we consider and characterize a wide class of black hole models that arise in 
the gravitational collapse of a massless scalar field, and which satisfy the regularity conditions 
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above. It is shown that for a large class of collapse scenarios the field collapses to a simultaneous 
singularity. 

We note that the second condition that we adopted above can be stated as 

M(0,i;) = ^ (51) 

for all V > 0, i.e. 1 > v > 0. At this point, it is useful to note that ^(i) is a free function 
and therefore can be chosen to be any function of t subject to its being regular. However, once 
the regularity conditions have been imposed, then two solutions (obeying the same regularity 
conditions), with two different functional forms of <j){t) need not be diffeomorphic to each other in 
general. 

We note that any r = constant curve is timelike, and the tangent vector of this curve is, 
= dx^/ds with components, 

r'' = (^,0,0,0) (52) 

The proper time along this curve is then given by r = / r^T^ds = J ds. Since ds'^ = e^'^dt^, in this 
case we get t = J e^dt, and we have, 

Mtf) 

T[v{tf),r] = / e^dt (53) 

Jv{ti) 

So t{v, r) is the proper time along any particular shell of comoving radius r to reach the value 
starting from the initial epoch v = 1. 

Subject to the above mentioned conditions, some general results can now be proved about the 
nature of singularities and black hole formation in gravitational collapse of scalar fields. We thus 
consider the classes of models where ?) < through out the evolution. Further, we take that v' >Q 
at r = 0, and at all other values of r. As we show here, this last condition would imply that 
this class of solutions would admit no non-simultaneous singularities as collapse endstate. Such 
a condition corresponds to the situation when the central shell at r = arrives at the spacetime 
singularity earlier in time as compared to other shells with greater values of r. This is related to 
avoiding the shell-crossing singularities R' = within the cloud which are not generally considered 
to be physically genuine. We shall show that, in such a case, the matter field then collapses to 
either a simultaneous singularity in a finite coordinate time, or that both ts{r) and the proper 
time T{v,r) must diverge along any timelike curve r = constant as the system evolves under the 
conditions we shall state below. 
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Proposition 1: If (j){t) is divergent at some instant ti, then there is a simuhaneous singularity 
at the time t = ti. 

Proof: The density p = ^e^^'^0(t)^. If is divergent, at t = ti, then there are two 
possibihties. The first possibihty is, the density remains finite and there is no singularity at t = ti 
in the case if we have e"^'^^'"'*^) = 0. However, this is not allowed because it means e'^'^^'^'^'^^ = oo, 
i.e. the metric component is blowing up at regular spacetime points. As this is not allowed by 
regularity conditions, the only other possibility that remains is the density must diverge at t = ti, 
so the singularity is simultaneous at the epoch t = ti. 



Therefore, it follows that the existence of a non-simultaneous singularity curve implies that the 
function (/>(t) must remain finite at all regular epochs in the spacetime. For a stiff fluid collapse, 
this means that if limn*^ ^ as limt — > ts, there must be a simultaneous singularity in the limit 
of t going to tg. 

Proposition 2: If tg{r) is not constant and if ti' > everywhere in the spacetime, then ts{r) 
must be divergent. 

Proof: First we note that from (I44p we get, 

v^{vM,,-3M-rM,r)^ 

""^^'^^^ wmm:. — ^''^ 

Therefore, 

We note that e'^^^'^^ when f 7^ by the regularity conditions. This implies that the other term 
in the denominator of the integrand must be finite at r = 0, otherwise ts(0) = 0, or the singularity 
will be present at the initial epoch itself, which violates the regularity of collapse from non-singular 
initial data. We note that M(0, v) = uiq/v^ for all values of v such that 1 > i; > as stated above, 
which is finite whenever v ^ 0. Therefore the quantity, 

li^J-|--^W--f^^-'-^^-^)'-ll,l,„X (56) 

must be finite. This finiteness condition and (j5ip together imply that M(r, v) must be of the form, 
M(r, f) = ^ + r"'g(r,v), where n > 2. This can be clearly seen by the direct substitution of this 
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form for M in ts(r) (see also Appendix A). Now we can write, 
liniX = lini^ -27^-1 +2r -t2^( n + 3 5 + 5,^ )) + ^^^^ ^/o T^I ^ i^'^) 

Now — 1) = fi{r), where /i(0) is a finite quantity in the Hmit of r — > because X cannot 

diverge. So we get, 

^ = 1 + ^^^^ 
When n > 2, we can see from above that liuir^oX = /i(0), and when n = 2, 

3 3 

hm^.^o-'^ = /i(0) + j7(o)9o{v) = /i(0) + ^go{v), where go{v) = {n + 3)g{0,v) + g,r {0,v). 

^From this, it follows therefore that in the case n > 2, 

lim lim vX(r, v) + Mir, v) = — ^ (59) 

When n = 2, 0[gQ{v)\ > 0(:;j7) which violates the continuity of Ts{r) and hence is ruled out (see 
Appendix B). Therefore we have, 



^1 

lim lim vX{r, v) + M(r, v) = ^ (60) 

r-*Ov—*0 V 



where ci is some constant, which is the result in general for n > 2. 

We now note that from (f^T]) . we have e'^^^''"^ = -^t^^- This can always be written as 

lim 6"^°'"^ = h{v) (61) 

The divergence of ts(0) if it is there can come from the range where v is close to zero. We denote 
that part of the integral by tsrf(O). So we can write. 



tsM = / ^-f-—rdv (62) 



Here C2 = rriQ or C2 = ci, and e is a small quantity. This can be written as. 



tsdiO) = I -^dv (63) 

x/C2 Jo Vh{v) 



Now, if the singularity is non-simultaneous then (j){t) remains finite always. Therefore, /3(f) is 
also finite. So we get the result that ts(0) is divergent. Since v' > 0, it follows that ts{r) is also 
divergent. This proves the required result. 

We note that if the singularity is simultaneous, then (/>(t) blows up when f (0, t) = by Prop. 
1. Then fz{v) must be divergent and then is(0) has to be finite. For a collapsing massless scalar 
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field, if the mass function Ai{r,v) and all the metric functions are at least near the central 
shell as they should be, and the singularity curve is non-simultaneous and an increasing function 
of time near the center, then the time taken for the central shell to reach the singularity diverges 
logarithmically, as we have shown above. 

As an illustrative example, we can consider the class of solutions within an e-ball around the 
central shell, where we can ignore the contributions of terms higher than even very close to 
f; — > 0. We note that this may not be possible for all classes of solutions, as in the Einstein 
equations there are terms of the form r"/t;™ {m,n > 0) and these may have non-zero limit in the 
vicinity of the singularity, that is w = 0. 

Since the mass function is near the center, without any loss of generality we can Taylor 
expand the function around r = as, 

M{r,v) = '^+M2ivy + ■■■ (64) 



We note that we get the zeroth order term by solving (j22p at r = 0. Also that the first order term 
should vanish follows directly from the no force condition at the center of the cloud. The function 
M2i'v) is an unknown function which solves (j22p near the center. Using (j64p we can write (|30p as 



v' = -^[v'^M2{v)],, 
bmo 



.A 



(65) 



Let us consider bo{r) = and /(r) = /(O) = 6mQ. Then putting the Taylor expanded forms of 
different functions in the (j2ip we get 



e-^ = l+('^^^^)r^ + 0(r«) (66) 



h{r,v) = {^-^^^^^]+0{r% (67) 



Also using ([S7|) we have 



3mo 

^From the above equations we get the behaviour of the function u near the center (since u' 
A^yR') as 

5v 



v'M2{v)Ur' + Oir'')+Cit), (68) 



12mo 

where C{t) is the constant of integration. Comparing with the densilty of the central shell we get 
C{t) = ln{y'^ /niQ). Therefore we can write, 



^3 

mo 



l + r^[v''M2{v)iy + 0{r^) 
\Zmo 



(69) 
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Finally using all the above equations and neglecting terms higher than the order r (since we are 
only concerned with an e-ball around the center), we get, 



\ + —r '■> = 

m-o V 3mo f 6mo 



+ ^;v' = -^[v'M2{v)l. (70) 



/,From (j32p we see that the integrability condition is trivially satisfied for r = 0. However, upto 
the first order correction we get the following condition, 



(71) 



The above is a complicated non-linear equation for the function A42(v). Let us consider that 
near the central singularity v{0, t) = 0, the function has a leading order power law behaviour 
M2iv) ~ v". Then using the above equation and comparing powers on both sides, we see that 
only consistant solution can be obtained with a = —3. In other words the behaviour of the 
function A42(w)near the central singularity is of the order of 1/f^. Using this in the singularity 
curve expression (|37|) we see that the time taken for the central shell to become singular diverges 
logarithmically. 

Returning to the class of collapse models under consideration, we shall now show that for this 
class non-simultaneous singularity cannot occur. First we define, 

ro{t)= f e-'^^'^Ut (72) 

This is the proper time along the r = shell. Now let us assume that To(tso) = tos is finite. ^From 
dZSl) we get, 

^ = e-^^O'*) (73) 
dt ^ ' 

We can then show the following result. 

Proposition 3 If > at tq = tqs, then for any > 0, Tr2{Tos) is divergent. 



Proof: One can make the following transformation which is allowed by the comoving coordi- 
nates, given by, : (t, r, 6, (p) — > x'^ = {To,r, 9, 0). In the new coordinate system, the line element 
is then written as. 

The new coordinates cover the manifold upto tq = tq^ only. Since, f ' > at tq = tq^, {r2,TQs) is a 
regular spacetime point. The proper time along the r2 = const, curve from tq = roi to tq = tqs is 
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given by, r,,(ro.) = Q ^M^)dro = ^ e'^^^-'^^^^'^^'^^^dro. 



Now v{0,t) = vo{t) is a function of r only. Therefore dt — i^(Kt)dvo — iifo^vn) "^^o- This implies, 



= T77^ rdvo (75) 



Using this, we get 



Here it is important to note that e'^l''^''"*'] = e'^V'^^'^i'^'^''"^)] jg finite and non-zero at vq = 0, because 
v{r2,0) > 0. 



^From the proof of Prop. 2, we have lim„g^o ^(0, wo) = —civq, where ci > 0. Since v goes to 
zero as ^ 0, a divergence may be present in r^j. So we consider that part of the integral where 
vo « 1, which is given by, Tr^{Tos)d = ^ /q ^^^^-^dvo, where e << 1. This gives Tr2{Tos)d = 
e [ 2.j3]Uo-o j.^e ^g^^p^ which is a divergent quantity. It follows that Tr2{Tos) is divergent. 

The proposition above proves that any non-simultaneous singularity cannot develop for this 
class of collapse models. That implies for this class only a spacelike singularity can form as collapse 
endstate. This characterizes a wide class of black hole formation models from a massless scalar 
field collapse, or for a stiff fluid collapse. 

We now also indicate here the possibility of existence of a class of non-singular solutions in 
evolving scalar fields, or stiff fluids scenario, by proving another proposition as below. 

Proposition 4: If v'{r, t) > b where 6 > 0, for ri < r < r2 for some ri, r2 > 0, and t G {U, oo), 
then we must have T{vi,r) > k for all A; > for all r, for some vi > 0. 



Proof: For r > r2, we have /J v'{r,t)dr = J^^ v'dr + J^^ v'dr + J^^ v'dr. 
This implies that, v'{r,t)dr > 5(r2 — ri) for r > r2. 

v'{r, t)dr = v{r, t) — f (0, t) 
Denoting 6(r2 — ri) by c, 
v{r, t) > v{Q, t) + c; for r >r2- 

".■ ■u(0, t) > for t G (tj, oo), therefore v (0, t) > c for f G {U, oo) for r > r2- 

For r > r2, the proper time elapsed during the comoving time interval (tj, oo), T{vi,r) = e^dt, 
for some vi > c. 

By the regularity condition, e'^^'"'") has a positive lower bound in the range v G (c, 1). 
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This implies that, r(fi,r) > k\/k > for r > r2 and for some vi > c. 

Since T{vi,r) is continuous in each t = constant hypersurface, and r2 can be any non-zero number 
however smah, there must exist some v > such that T{v,r) > k\fk > for r < r2- So T{v,r) > 
k'\/k > for Vr for some v > 0. This proves the result. 

Thus all the classes of solutions that satisfy the conditions assumed as above would be singularity 
free. If there are no such solutions, then that would indicate the possibility of bouncing models 
within the framework above. We note that since singularities would not form in this class of models, 
there would not be any formation of trapped surfaces also. We discuss this in some detail below. 

For spherically symmetric spacetimes, the equation of the apparent horizon is given by 

g>'''R,^R„ = (77) 

In this case, this implies G — H = 0. Using ([8j), one can rewrite the condition for formation of 
apparent horizon as F = i? or as 

M = ^ (78) 

Proposition 5: For the non-singular class of models discussed in Proposition 3, the apparent 
horizon does not form in any finite coordinate time. 

Proof: To prove this, we first prove another result. 

Lemma A: We must have mod v < 5, ^6 > 0, in a, total time interval, which is infinite. 
Proof: We know, ts{r) = oo, Vr. Now we have, 

rts{r) 

v{r,ts{r))-v{r,ti) = / vdt. (79) 
Jti 

This can be written as 1 = f^{—v)dt = J^*^ {—v)dt + f^^ {—v)dt + = Ii + 12 + where (ti —ti), 

(^2 — ti), ....are all finite. 

This implies that, given any 61 > 0, there exist infinite number of integrals /„s such that, In < 61. 
In each of these integrals, mod v < S, \/d > 0; i.e. mod(change of v)< any arbitrarily small 
amount. Since there are infinite number of such integrals, mod 1) < 5, V(5 > 0, in a total time 
interval that is infinite. (Proved) 

Since mod v < 6, ^6 > 0, for a total time interval which is infinite, there must exist a time 
instant ti at which this holds, such that, ti < tgo and ti > I for any given / > 0. 
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■-■ is{r) > tsOi "^f-, so for some v > 0, v{r,ti) > for r > 0. By the regularity condition, at i = ti, 
^^{■r,ti_) y Q From Einstein equations we have ^/vv = —e'^[^{G — 1) + M]^ . 
Therefore at t = ti, 

[^{G - 1) + M] <5,y5> 0. (80) 

By the energy condition, p = -M^ > 0. This imphes M,„ < 0. From i^, G > 0. 

Let, v{ti,r) = Vf{r). Then the inequahty can be written as M < 6 + {1 — G)^, V(5 > 0. But we 

know, (1 — G)^ < This imphes 

M{r,vf)<^-1. (81) 
Let Ve{r) be the value of v{r) at any time earlier than ti. 

^P- < ^(Since v < 0). Also M{vf,r) > M{ve,r). These imply M{ve,r) < ^ < ^ or 

Therefore the equation of apparent horizon (j78p is not satisfied for r ^ and when ti < t < ti. So 
there is no apparent horizon and consequently no trapped surface formation in that time range. 
Since there always exists some ti such that ti > I for any given /, for r / 0, trapped surface does 
not form in any finite time. 

For r = 0, M = Putting this in ([78]) . mor^ = v^. This is satisfied when v = 0, i.e. at 
t = tso- Since tso = oo, apparent horizon does not form at the r = shell at any finite time. So for 
such a solution of the Einstein equations, the apparent horizon or trapped surfaces do not form at 
any finite coordinate time. (Proved) 

In Proposition 4, it has already been shown that an infinite coordinate time interval correspond 
to infinite proper time interval along any r = constant world line. This implies that for this class 
of models, no apparent horizon or trapped surface would form at any finite proper time. 

C. Models collapsing to form non-simultaneous singularity 

Now we discuss briefiy the other classes of models for which the singularity can be non- 
simultaneous. Toward such a purpose, we must clearly relax the condition v' > 0, thus allowing v' 
to be either positive or negative. All the other regularity conditions still remain the same. 

In the earlier class of models, the divergence of ts{r) was responsible for the fact that there 
could only be simultaneous singularity. The divergence of M as — > caused the divergence of 
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ts{r). So one must look for such a form of M which is different in this respect. One such generic 
form of M is given by 

where n > 2, gi{r, 0) blows up and 52 0) is finite. In this case M(r, 0) would remain finite except 
at r = 0. 

Here ts{0) will be divergent as in the previous case, because M{0,v) = But in this case, 
ts{r) would be finite for any r > 0. To show this we first show that there is no divergence in ts{r) 
when evaluating the integral near v = 0. We know that lim^^Q M{r,v) = r^g2{r,v), where 52 is a 
regular function. Also, limj;_^o -^(^; ^^)) = f''^g2{f,v),v and 

lim.- = -'^^^^'°)+:^^^:'°)'^ (83) 
v^o 2rg2[r,v),v |„=o 

^Prom ()44p , it is seen that limt,_,o G ~ and from (j4ip , lim^,.^ Putting them all in 

([27ll . we have limt,^o v ~ —\/v- This means that there would be no divergence in the ts{r) integral 
for any value r > coming from the part near v = 0. So ts{r) is finite for r > 0. We note that 
since ts{0) = 00 and ts{r) for r > is finite, f ' > is not satisfied. The proper time for the r = 
shell to reach the singularity is, however, finite as we can show. The proper time, r^o = e'^dt 
can be written in terms of an integral of v also. Then r^o = Ji ^dv. If we consider the part of 
the integral where v « 1, it is given by j^v'^dv, where e << 1. This goes to zero. So the proper 
time TgQ is finite. Thus for the above generic form of M, a non-simultaneous singularity may occur. 
Such a singularity may also be spacelike or timelike in this case. 



V. APPENDIX A: GENERAL FORM OF M 



Since M(0,t;) = M can in general be written in the form, 

M{r,v) = '^ + f2{r)g{r,v) (84) 

where g{o,v) + and h{^) = 0. Then M„ = + h{r)g,, and M„. = f^{r)g{r,v) + f2{r)g,r. 

iFVom (I22D, 

, v'^{3f2{r)g + rf^{r)g + rf2{r)g,r +vf2{r)g,v ) 



2r(3mo - f2{r)v^g,v ) 

This gives 



(85) 



T,/ _ , _ ,,ri , ^'(3/2(r)g + r/^(r)g + r/2(r)g,, +vf^{r)g,,) 

^ 2(3mo-/2(r)A,.) ^ ^ ^ 
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^Prom (j44p . we have 

= + v^{3f2{r)g + r/2(r)5r + rf2{r)g,r +vf2{r)g,^ ) + 

(3/2(r)g + r/^(r)g + r/2(r)ff,, +z;/2(r)g,. )^ 
4 (3mo-/2(r)A,.) ^ ^ ^ 

For regular epochs (i.e. v ^ 0),M and its derivatives are bounded. In the integrand of ts{r), we 
have in the denominator, a factor [^{G — 1) + M]^ , which must be finite when r goes to zero. 
Since M(0, v) is finite, ^{G — 1) must be finite. This means limr — > 0^(G — 1) is finite. We can 
write the quantity ^(G — 1) as, 

1 Qttiq 1 '2v'^ 

- 1) = ( ^2j2(^) -^) + ;272(;:y (3/2(^)5 + rf^{r)g + rf2{r)g,r +vf2{r)g,, ) + 



(3/2(r)<7 + rf^{r)g + r/2(r)5„ +T;/2(r)<7,. )2 



2r2/2(r) (3mo-/2(rK5„) 
Now the first term in the bracket and the coefficients of and have to be zero separately. This 
gives, 

1 + r^ji(rj 

The coefficient of is g^(/i(r) + ^)(3/2(r)s- + r/^(r)g + r/2(r)5r,r +t!/2(r)5(,„ ). This implies that 
limr 0^^^^ and limr 0^-^^^|^ are bounded. This implies that as r goes to zero, f2{r) goes to 
zero atleast as fast as r^. 



VI. APPENDIX B: RULING OUT SOLUTIONS FOR WHICH 0[goiv)] > 0(i) 

4 

In this case, limt,_^o[^(G - 1) + M]r=o = ^go{v). Let, liuiy^o go{v) = where 6 > and 
C3 is some constant. Using this and ([6T]) . we get, tsd(O) = — rdv. This can be written 

as tsrf(O) = y^^b^lo- This is vanishingly small. So ts(0) is finite. 

However, for r > 0, lim^^o[^(G- 1) + M]^=o = ^ and lim^,_,o e'^^'^'^^ = t;^+i/3(r, w). Putting 
them together, we get, tsd{r) = ■;^jj^y^dv. This implies that ts(r) is divergent when r > 0. 

Now there are the following possibilities. Firstly, the singularity can be non-simultaneous, i.e, 
t;' > at the epoch when some r > shell reaches the singularity. In this case, since tg{r) = 00 
for r > 0, arguing in the same way as in Prop. 3; it can be shown that Tg{r) will be divergent for 
some value of r. This means that the singularity curve Ts{r) is discontinuous. So this possibility is 
ruled out. 

There is also another possibility that all the shells r > hit the singularity at the same time, 
i.e. v' = Q for r > at that epoch. Since = at the epoch of singularity of some r > shell, 
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all r > shells reach the singularity at the same epoch. The proper time taken along a shell with 
a very small value of r would be close to the proper time taken along the r = shell to reach 
the singularity. If the Ts(r) is continuous, then it follows that at tq = r^o, all the shells reach the 
singularity. But this contradicts the fact that ts{0) < ts{r) for r > 0. So this second possibility is 
ruled out. 
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